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dimensional bulk. At low energies, the effect of the bulk on the brane evolution 
can be described through a mirage component, termed generalized dark radiation, 
in the effective four- dimensional field equations. Using the covariant formalism, we 
derive the exact form of these equations. We also derive an effective conservation 
equation involving the brane matter and the generalized dark radiation. At low 
energies the coupled brane-bulk system has a purely four-dimensional description. 
The applications of the formalism include generalizations of the Starobinsky model 
and the Dvali-Gabadadze-Porrati cosmology. 
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1. Introduction 

The novel characteristics of braneworld cosmological models in the context of the 
Randall- Sundrum (RS) scenario || can be understood in terms of two basic ingredi- 
ents: the effective compactification of gravity near a brane within an Anti-deSitter 
(AdS) background, and the appearance of a dark (or Weyl or mirage) radiation term 
in the effective equation describing the cosmological evolution [Q]. This term results 
from the presence of a black hole in the AdS bulk spacetime, whose gravitational 
field affects the brane motion y[. 

Starting from this setup, several generalizations have been considered, which 
include modifications of either the 5D bulk spacetime or the brane gravity sector. For 
example, the leading higher- curvature bulk correction in the context of string theory 
takes the form of a Gauss-Bonnet (GB) term. The inclusion of such a term results in 
an extension of the conventional braneworld scenario with unusual features ||, such 
as the existence of a self-accelerating branch of the cosmological expansion. However, 
this branch suffers from ghost-like instabilities, which make its physical interpretation 
problematic ||. Radiative corrections can also modify the 4D brane action, leading 
to the Dvali-Gabadadze-Porrati (DGP) model || and its generalization for non-zero 
brane tension and bulk cosmological constant @. The combination of the DGB 
model with a bulk GB term has also been considered ||. The cosmology of the DGP 



model has similar properties and pathologies [|10j as the braneworld cosmology 
with a bulk GB term. 

Another direction in the generalization of the braneworld scenario is to include 
matter in the bulk action and analyze the effect of its gravitational field on the 



brane dynamics, as well as its direct interaction with the brane matter [IT]. Finding 
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exact and general enough solutions for the bulk geometry is a difficult task. In the 
original studies the influence of the bulk matter on the brane cosmological evolution 



was tractable only in some simple cases [12, 13 1. This difficulty was overcome by 



employing covariant techniques in order to describe the solution of the braneworld 
equations. The covariant approach was used in [TJ|, [TJ|, |16[ , where it was found that 
the solution depends on the so called "comoving mass" Ai of the bulk fluid. This is a 
function of the brane scale factor R and the proper time r on the brane. In the case 
of an AdS-Schwarzschild bulk it is constant, equal to the black- hole mass. The form 
of Ai (R, r) can be determined only within a specific model of the bulk dynamics. 
Such a model may involve several bulk fields that possibly interact with the brane, 
or may employ a description in terms of a bulk cosmological fluid with a certain 
equation of state. The exact form of Ai(R, r) is necessary for a detailed discussion 
of the cosmological evolution of the brane. Nevertheless, the general properties of 
the evolution, such as the occurrence of accelerated expansion, or the rate of energy 
exchange between the bulk and the brane, can be determined by the structure of 
the resulting covariant equations i.e. the Friedmann and Raychaudhuri equations 
[|l7jl . The "comoving mass" appears in the Friedmann equation within a term ~ 
A4(R,t)/R 4 , characterized as generalized dark radiation. 

Recently, there have been proposals to consider f(TZ) generalizations of the Ein- 
stein term in the bulk or brane action. The inclusion of such a term in the bulk sector 
IS, 1J| is possible. However, the covariant definition of the junction conditions [20 



on the brane becomes ambiguous because of quadratic 5-function contributions to 



the field equations \TE\. Only the Gauss-Bonnet combination of higher derivative 
terms results in well defined junction conditions [||. The inclusion of higher deriva- 
tive terms on the brane, represented by a smooth function f(TZ) of the brane scalar 
curvature, does not lead to such difficulties |2l], Moreover, it provides a direct 
generalisation of the DGP model. The purpose of the present paper is to give the 
set of equations describing the most general formulation of this problem: brane evo- 
lution with higher derivative contributions resulting from a brane f{TZ) term, in the 
presence of arbitrary matter components on the brane and in the bulk. 

In section 2 we present some general and covariant results regarding the geomet- 
rical and dynamical structure of the bulk-brane setup. We also give the exact form 
of the induced Friedmann and Raychaudhuri equations. As expected, apart from the 
higher derivative corrections, these equations contain the generalized dark radiation 
term. We also show that there is conservation of energy between the brane matter 



and the generalized dark radiation, generalizing in this way the result of [|14[ |15 



In section 3 we solve the field equations in a specific coordinate system by assum- 
ing a static and anisotropic fluid in the bulk. We reproduce the Friedmann and 
conservation equations. In section 4 we show that, at low energies, the brane cos- 
mological evolution is that of a four-dimensional Universe. The evolution possesses 
a self-accelerating branch exhibiting three new properties: a) an effective cosmologi- 



-2- 



cal constant, b) an effective negative energy density associated with the generalized 
dark radiation, and c) effective f(TZ) corrections. Finally in section 5 we draw our 
conclusions. 

Throughout this paper the following conventions are used: The 5D bulk space- 
time manifold is endowed with a Lorentzian metric of signature (— , +, +, +, +), bulk 
5D indices are denoted by capital latin letters A, B, ... = 0, 1, 2, 4 and greek letters 
denote brane indices a, (3, ... = 0, 1, 2, 3. 



2. General covariant results 

The effective action we consider has the form 

S = J d\^~g (A + M 3 TZ + CZt) + f d'x^g [-V + CZl e + r c M 3 f{K) 

(2.1) 

In the first integral, —A is the bulk cosmological constant (we assume A > 0). In the 
second integral, V is the brane tension, g a p the induced 4D metric on the brane, g its 
determinant, TZ the corresponding curvature scalar, and r c the characteristic length 
scale of induced gravity. The matter contributions are arbitrary, and the effective 
action incorporates possible quantum corrections in this sector. We assume, however, 
that the corresponding energy-momentum tensor is consistent with the underlying 
geometry. 

The Einstein field equations (EFE) take the form 

° AB = 2lP ^ + ^ 9AB) ' (2 - 2) 
with the total effective energy-momentum (EM) tensor T A °^ given by 

T a *=T ab + 5( V )t ab . (2.3) 

The term T A b is the bulk matter contribution, while r AB is the contribution from 
the brane located at rj(x A ) = 0. The presence of the higher derivative gravity term 
in the brane action, represented by the function f(li), implies that the tensor r a p 
takes the following form [^2] 



r a p = f aP - V~g aP - 2r c M 3 S Q/3 . (2.4) 

Here T a p is the brane EM tensor and the function f(7Z) results in a contribution to 
the tensor r a p proportional to defined as p3| 



£ Q/3 = f'(K)Kafi - \f{K) ~9ap - VaVpfXK) + g a ^ 2 f'(il), (2.5) 



where f'(TZ) = df/dTZ and the covariant derivatives are constructed with g, 



a/3- 
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In the analysis of cosmological setups, it is natural to restrict our formalism 
to geometrical backgrounds with a high degree of symmetry. We make the sim- 
plest symmetry assumption and derive the general form of the solution of the above 
equations in the case of a Friedmann-Robertson- Walker (FRW) brane. As a result, 
there exist 3D hypersurfaces T>, invariant under a six-dimensional group of isome- 
tries. This implies that the surfaces T> have constant curvature, parametrized by the 
constant k c = 0, ±1. It should be noted that the assumption of maximal symmetry 
of T> implies the existence of a preferred spacelike direction that represents the local 
axis of symmetry, with respect to which all the geometrical, kinematical and dy- 
namical quantities are invariant. This direction can be chosen in two different ways 



14], which are observer-depended: either w.r.t. the (prolongated) brane observers 
u A , for which the spatial direction is represented by the vector field n A {n A riA = 1, 
u A riA = 0); or w.r.t. the (comoving) bulk observers u A for which the spacelike vector 
field is e A [e A tA = 1, u A eA = 0). The two approaches yield identical results for the 
cosmological evolution. Because of the presence of the term in the total brane 
EM tensor, we find it convenient in the present paper to apply the methodology 
of Hl6l , which is based on the comoving brane observers. In order to simplify the 
analysis, we make the usual assumption of a Z 2 symmetry around the brane. 

The structure of the bulk-brane setup can be studied using elements from the 



theory of spacelike congruences. We introduce the projection tensor Jl6| 

IIab = 9ab + u A u B ~ n A n B = h A B ~ n A n B (2.6) 

n/ = 3, u A u B c = u A , u A n B = u A u B = o, (2.7) 

which is identified with the associated metric of the 3D manifold T> (the screen space), 
normal to the pair {u A ,n A } at any spacetime event. The first covariant derivatives 
of the spacelike vector field n A are decomposed according to 

n A - B = - n A u B + n A n B = K AB + n' A n B . (2.8) 

Here 

■& = n A -, B ^ AB = n A A + n A u A (2.9) 
is the rate of the surface expansion V, and we have used the notation 

P'a... = Pa..,lu l , P a ... = Pa..,lu l (2.10) 

for the directional derivative along the vector fields n A ,u A of any scalar or tensorial 
quantity. In addition Kab = (^/3)n^ B — uaUb corresponds to the extrinsic curvature 
of the timelike hypersurfaces normal to n A . 

The definition of the overall expansion of the spacelike congruence implies 



ii 



■d = U AB K AB = n\ A = 3 T , (2.11) 
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where I can be interpeted as the average length scale of V. For example, in the 
spherically symmetric case k = 1 it represents the radius of the spheres T>. The 
temporal change of I is controlled by the expansion rate of the timelike congruence 
as measured in the screen space V, namely 

U AB u A , B = u\ A = 3~. (2.12) 

Therefore, the length i completely determines the volume of T>. For on brane con- 
siderations, it is identified with the scale factor of the brane Universe. Employing 
the Gauss-Codacci equations for the spacelike surfaces T>, it can be shown that [ffBl 





Yi=-\{ E + \ G ^ AB - 2G -)-\l \ +(t) (2-13) 



7 I + I 7 1 =\{ £ + \ G *B9 AB - 2G±) - jn A u B - \c AB n A n B . (2.14) 

The generalized dark radiation term £ = C A cBDU A n c u B n D = C A cBDU A e c u B e D is 
the spatial eigenvalue of the electric part of the 5D Weyl tensor, and 3G± = H AB G AB . 
Equation (|2.13|) shows how the scalar curvature of the 3D space T> is affected by the 
kinematics and the dynamical (when the 5D EFE are employed) content of the 
spacetime. At the location of the brane, equations (|2.13|) and Q2.14|) correspond to 



the effective Friedmann and Raychaudhuri equations of the FRW brane. 

The matter content of the bulk is described by the EM tensor T AB , which can 
be written in the usual way with respect to the bulk observers u A : 



where 



T AB = pu A u B + ph AB + 2q( A u B ) + n AB , (2.15) 



p = T AB u A u B , p = -T AB h AB , q A = -h c A T CD u D , (2.16) 
vtab = h c A h B T CD - hi CD T CD h AB . (2.17) 



The first term in the rhs of equations ( |2.13j ) and ( |2.14| ) can be determined through 



the full 5D EFE ([Oft and the bulk EM tensor ( pT5|) . The result is ||, |15 

k_ M A_ _ 

£2 ~ 6MW 12M 3 




(2.18) 




I'. ~a M A 1 

'T AU ~ 6MW " 12M3" 6M^"' (2,19) 
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where pn = Tab^ A ^ b is the pressure along the direction of the preferred spacelike 
vector field n A , and M. is the "comoving mass" of the bulk fluid, satisfying 

(M - Mo)' = 2n 2 p£ 3 £'. (2.20) 

Only in the spherically symmetric case (k = 1) the "comoving mass" M. has the 
usual physical interpretation as the effective gravitational mass contained within a 
sphere with radius £. However, we shall refer to M. as the "comoving mass" for all 
geometries of the hypersurfaces T>. The integration constant A4q in equation ( |2.2C| ) 
can be interpreted as the mass of a black hole at £q = 0. 

The covariant form of the junction conditions for braneworld models has been 
derived in [B3J] and involves the projected (perpendicular to n A ) part of the extrinsic 



curvature 

K *P = -4^3 ( r -/3 - \ T 9*pj ■ ( 2 - 21 ) 

Equation ( p.21| ) can be used in order to determine the discontinuous quantities £'/£ 
and tiau a by projecting along Ii a ^ and u a u^ — 2Il Q ^/3 respectively. In addition, 
assuming a perfect fluid matter configuration on the brane, the brane EM tensor can 
be written in terms of the energy density p and the isotropic pressure p as 

f a p = pu a U/3 + ph aP . (2.22) 

At the location of the brane, we deduce that 

ws{ v+ ~p- 2r < M3 [ 3 f ( H2 + fa" -f)+ 3fl£r] } 



12M 3 

n A u A = {3p + 2p-V + 2r c M a 



9/'| + 3 (/')'• -nf + f 



(2.23) 
(2.24) 



We recall that for on-brane considerations, £ = R, with R the scale factor of the 
brane. The directional derivative denoted by a dot (defined in equation |2.10|) becomes 
a derivative with respect to the proper time on the brane. As a result, H = R/R is 
the rate of the cosmological expansion. Also, /, /', /" are considered functions of 
1Z = 6(H + 2H 2 + k c /R 2 ). The sign of the directional derivative £' determines the 
way the brane is embedded in the bulk spacetime. For r c = 0, a negative sign for £' 
results in localization of the low energy gravitons near the brane. It is apparent from 
equation ( |2.23|) that £' < corresponds to a brane with positive tension, a property 
that guarantees stability under small perturbations. The self-accelerating branch of 
the DGP model |g || at late times has V = 0, f{K) = K, p, a' 1 -»■ and H 2 ~ l/r c . 
It is then apparent from ( |2.23[ ) that £' > in this case. This branch is known to have 
ghost-like instabilities 
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Using equations ( 2.18 ) and ( p. 23 ), the effective Friedmann equation becomes 



H 2 + ^; + 



A M(R,r) 



R 2 12M 3 Qn 2 M 3 R 4 



144M 6 



V 



2r r M s 



3HTZf" 



(2.25) 



This equation represents the most general form of the brane cosmological evolution 
in the presence of bulk and brane matter and higher derivatives terms in the brane 
action. The bulk affects the brane evolution through the "comoving mass" M, which 
is a function of the scale factor R, but may also have an explicit dependence on the 
proper time on the brane. The effective Friedmann equation ( |2.25| ) results from 
squaring £'. It includes both the normal and self-accelerating branches of the DGP 
model. 

Throughout this paper we assume the standard tuning to zero of the effective 
cosmological constant in the Randall-Sundrum scenario This is achieved if the 
bulk cosmological constant —A and the brane tension V are related through A = 
V 2 /(12M 3 ). We also define the energy scale k = V/(12M 3 ) = [A/(12M 3 )] 1/2 . The 
term 

M(R,t) 



h = (2 ' 26) 

is the effective energy density of the generalized dark radiation [I7| and reflects the 
influence of the bulk matter on the brane evolution. 

The Raychaudhuri equation ( |2.19|) can be used in order to derive the effective 
conservation equation between the brane and the bulk. Taking the derivative of the 
Friedmann equation ( 2.25 ) along u a , combining the result with (|2.19 ), and using 
( |2.18|) , fl2.23| )-( [2.24|) , we find the remarkably simple effective conservation equation 



p + 3H(p + p) 



-e 2 k 



Pd + H ( 4pd + -j-P\\ 



k k N ~ 1/2 

H2 + w + k2 ~6W^ 



(2.27) 



The two possible values e 2 = 1 or — 1 correspond to the two embeddings of the brane 
in the bulk spacetime. For r c = 0, these require a brane tension V > or V < 0, 
respectively. 



3. A case study: Anisotropic bulk fluid 

In the previous section we derived the most general form of the Friedmann and 
Raychaudhuri equations for an arbitrary bulk matter configuration. However, it is 
instructive to analyze a class of cases in which the problem is tractable in specific 
coordinate systems. We assume that for a certain observer the bulk content can be 
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described as a static fluid. This assumption allows the possibility of an arbitrary 
number of fields and relies only on the existence of an observer comoving with the 
bulk matter. Clearly, important physical situations, such as those that involve the 
propagation of electromagnetic or gravitational radiation, are excluded by our as- 
sumption. However, many interesting backgrounds, including generalized black-hole 
ones, are allowed. 

In order for the embedding of a cosmological 3-brane to be possible, the spatial 
part of the metric must include a 3-space of constant curvature. The resulting metric 
can be cast in the form 



The lhs of the EFE 



ds 2 = -n 2 (r)dt 2 + r 2 dn 2 + b 2 (r)dr 2 . 
take the form 



3 1 



V 



G° = ^(---)-^ 



3 k r 









-I 




H2i) 


l 1 - 








r 


\r 


n J 


b 


^ n 


r / 



n 



- 



3 1 



b 2 r \r 



n 



G 4-7^- - + — - -o-, 



n 



3 k r 



(3.1) 

(3.2) 
(3.3) 
(3.4) 



where here the prime denotes a derivative with respect to r. 

The general form of the bulk EM tensor consistent with the above geometric 
setup is 

T A B = diag(-p,p,p,p,p) , (3.5) 

with the two pressures p, p not equal unless the bulk matter can be interpreted as a 
perfect fluid. The 00 component of ( |2.2|) gives 



2k c r + 



3M 3 



r 3 (A-p), 



(3.6) 



whereas the combination of the 00 and 44 components results in 
The conservation of the bulk EM tensor can be written in the form 



3(p-p) 



1 



p + p r(p + p) 

Because of the Bianchi identities, the set 
Integrating (|3.6| ) we find 

,2 



6M 3 



(p + A)rb 2 + 



Kb 2 



V 2 Ar 4 
V = cT + 12M^ 



(3.7) 
(3.8) 

1|) completely describes the solution. 
M{r) 



6vr 2 M 3: 



(3.9) 
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where Ai(r) satisfies 



dM 

dr 



2it 2 r 3 p 



(3.10) 



and corresponds to the "comoving mass" of the bulk fluid. 

In order to analyze the cosmological evolution of the brane, we employ the Gaus- 
sian normal coordinate system in which the metric takes the form 



ds 2 = —m 2 (r, i])dr 2 + a 2 (r, r])dVL 2 k + dr] 2 , 
with m(r,T] = 0) = 1. Through an appropriate coordinate transformation 



t = t(r,r)), 



r(r,r]) 



(3.11) 



(3.12) 



the metric (|3.1|) can be written in the form of equation ( |3.11|) . We define R(r) = 
a(r, r] = 0). In the system of coordinates (t, r) of equation (|3.1| ) the brane is moving, 
as it is located at r = R(r). Hence [|12| 



dt_ 

m 

dr] 
da 

da 
dr] 



n(R) 

m 

n(R) 



b 2 {R)R 2 



1 



1/2 



-e 2 



R 



R 



-£■2 



b(R) 



b 2 {R)R 2 + 1 



1/2 



(3.13) 

(3.14) 

(3.15) 
(3.16) 



where the dot denotes a derivative with respect to proper time and e 2 = ±1. The 
rj- derivatives are evaluated for rj = + . The value of e 2 determines the way the 
brane is embedded in the bulk space. As we have mentioned earlier, we impose a Z 2 - 
symmetry around the brane. We consider a matter configuration that is consistent 
with the solution of the EFE in an infinite bulk before the brane embedding. When 
the brane is included, only the solution in half of the space and its mirror image are 
employed. The value of e 2 determines which half-space is used. A negative sign in the 
rhs of equation ( 3.161) means that r decreases away from the brane. In the absence 
of induced gravity, the brane has positive tension. The configuration is stable under 
small perturbations and the massless graviton is localized near the brane. 

The bulk EM tensor at the location of the brane in the coordinate system (r, rj) 
is (no summation over repeated indices) 



T 00 = p(R) + [p(R) + P(R)} b 2 {R)R 2 
T tl = R 2 p(R) 

T 44 = p(R) + [p(R)+p{R)) b 2 {R)R 2 



T 04 = e 2 b{R)R b 2 {R)R 



1/2 



\p(R)+p(R)}. 



(3.17) 
(3.18) 
(3.19) 

(3.20) 
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The sign of T04 indicates whether a brane observer detects inflow or outflow of energy. 
This sign is determined by the value of e 2 and the rate of expansion (or contraction) 
R/R of the scale factor on the brane. 

The lhs of (|2.2|) near the brane (77 — > ) take the form (no summation over 
repeated indices) 

3a' 2 

-,2 



G { 



G 1 



G 



G { 



a 



[a 



3k c 

IF 



3d 2 3a" 

+ 

a^ a 

2m'a' a 2 

+ ^ + 

a* a a z 

12am' (a') 3 3 (a n2 



a- 

3m'd 3d' 

+ — , 

a a 



2a" „ 
h m — 


2d 


k c 






a 


a 


~F 


3m'a' 


3d 2 


3a 


- + 






a 


"a 5 " ~ 


a 



3k r 



(3.21) 

(3.22) 

(3.23) 
(3.24) 



with a prime now denoting a derivative with respect to 77. 

We consider a brane Universe containing a perfect fluid with an EM tensor 

m 2 (r, 77) „ 



<5(r/)a (r, ?])diag 



a 2 (r,r/) 



■p,p,p,p,0 



(3.25) 



Integrating the 00 component of ( |2.2| ) on a small 77 interval around the brane and 
using 



a' 



+ 
a 



we obtain 
1 



12M 3 



V 



P 



2r r M 3 



3f 



a 1 



\(n.r-.f -wnr 



(3.26) 

where H = R/R and /, /', /" are considered functions of TZ = 6(H + 2H 2 + k c /R 2 ). 
From ( |3.26|) and ( |3.16| ) it is straightforward to derive the effective Friedmann equation 

In the low-energy limit (p, H, a^ 1 — > 0) our choice of sign for a' + in (|3.16| ) must 
be consistent with the rhs of ( p.26| ) . For example, for r c = a negative sign (62 = 1) 
for a', is consistent with the negative sign in the rhs of ( |3.26| ) only if V > 0. The 
self-accelerating branch of the DGP model || [J at late times has V = 0, f(7Z) = TZ, 
p,a~ x — > and H 2 ~ l/r c . It is then apparent from (|3.26| ) that a' + > (e 2 = — 1). 
The effective Friedmann equation ( p.25|) results from squaring a' + . As a result, it 
includes both the normal and self-accelerating branches of the DGP model. 

The effective Friedmann equation can be complemented by a conservation equa- 
tion for the energy densities that drive the expansion. In order to obtain this equation 
we integrate the 00 and ii components of ( |2.2j ) on a small 77 interval around the brane. 
This gives (no summation is implied) 

1 n 



a + 
a 

m'. 



12M 3 
1 



4M 3 



T 



6M 3 



(3.27) 
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at the location of the brane. From equation ([3.24 ) we now find 



4M 3 



[r° + 3H(r° - r\ 



(3.2* 



at the brane. The modified energy momentum tensor of the brane, given by equation 



(234]), contains a piece H a p that arises from the induced gravity term in the action 
(2.1). We have verified that H a p is conserved. Therefore, the rhs of equation ( j3.2S|) 
receives contributions only from the density and pressure of the brane matter, so 
that 



G c 



AM 3 



(3.29) 



The bulk energy momentum at the location of the brane, as measured by an ob- 
server comoving with the brane, is given by equations ( |3.17| )-( |3.20| ). Solving equation 
( P-19|) for p(R) and substituting in equation ( |3.20|) gives 



T° 4 = -e 2 ^[T u + p(R)) 



1 



+ 



R 2 



R 2 b 2 (R) R 2 



-1/2 



(3.30) 



The energy density p in the above expression corresponds to the bulk energy density 
as measured by a bulk observer. It is more convenient to express it terms of the 
energy density of dark radiation ( |2.26|) . Making use of equation ( |3.10|) we find 



P(R) = ^ (P d + 4Hpd) ■ 



(3.31) 



The component T44 determines the bulk pressure as measured by a comoving brane 
observer. It is equal to the quantity p\\, covariantly defined below equation ( |2.19| ). 
We define an effective pressure term for the dark radiation through the relation 



Pd 



T + 3k Tu - 



(3.32) 



Making use of these definitions, the 04 of the EFE at the location of the brane, and 



the solution of equation (|3.9| ) for b(R), we find the conservation equation 



-[p- + 3H(p + p)} = 



M Y 1/2 r 



(3.33) 



4. Effective equations at low energies 



The main equations that have resulted from our study are the effective Friedmann 
equation ( |2.25| ) and the conservation equation ( |2.27| ). The Friedmann equation can 
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also be written as 



-e 2 



1 + kTc f + 2kr c ^ + 2kr c ^ - 2(kr c ) 2 ^ 



1/2 



• (4-1) 



The term 



A = -6r c M 3 



(/' " 1) ( H 2 + ^) - I (nf - /) + HTZf" 



(4.2) 



incorporates the effects of the deviation of f(TZ) from the standard Einstein form 

f = n. 

For p, pd, A <tiV, keeping only the linear terms, we obtain 



H 2 + — + 1 ~ 62 + kTc 
R 2 1 + kr c 

2(1 -e 2 )(l + kr c ) l-e 2 + kr c p 



(/' - 1) (h 2 + ^)-\ (£/' - f) + ff£r 



+ 



£2 



fcr c (l + A;r c ) 6 (M 3 /k) l + kr c 6 (M 3 /k) ' 



(4.3) 



with TZ = Q(H + 2H 2 + k c /R 2 ). 
For e 2 = 1 we have 



H 2 + 



6M 2 Y 



Ay 



+ 



kr r 



R 2 1 + kr c 



Ay 



i? 2 



(/' - 1) [H 2 + -£ - - [TZf - f + HTZf" 



6 



(4.4) 



with M|j = M 3 (r c + l//c). For A;r c 3> 1, the lhs takes the exact form it would have 
for four-dimensional f{TZ) theories. The sources, including the energy density of the 
generalized dark radiation contribute linearly in the rhs of the same equation. For 
kr c <C 1, the lhs of the Friedmann equation takes the conventional form, as expected 
for the low-energy limit of RS cosmology. 
For e 2 = — 1 and kr c > 1 we obtain 



f'[H 2 + ^\ -l[Kf-f)+ HTZf" 



R 2 



6 



with Mpj = M 3 r c . This is the Friedmann equation in f(TZ) gravity with two new 
features: a) An effective cosmological constant appears, despite the fine tuning of the 
bulk cosmological constant and the brane tension, b) The other striking feature is 
the negative sign of the contribution proportional to the energy density of the dark 
radiation. 
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For 62 = 



— 1 and kr c < 1 we have 



4 1 

^7 + 



kr, 



tt2 i 



-^[nf'-f)+HTZf" 



p 



-pd 



(4.6) 



with iWfj 



M 3 r c /2. 

It is obvious from the above that the brane cosmological expansion in the branch 
with €2 = — 1 has novel properties arising from: a) an effective cosmological constant, 
b) an effective negative energy density associated with the generalized dark radiation, 
and c) effective f{TZ) terms, that persist even in the limit r c — > 0. The second feature 
is not a consequence of the violation of the weak energy condition, as the energy 
density is assumed positive both in the bulk and on the brane. The first two features 
also appear in the DGP model |J, characterized by A = V = and / = 1Z. 

In the case of an AdS-Schwarzschild bulk we have pd ~ R~ 4 . At late times the 
contribution from the dark radiation is subleading to the contribution from the brane 
matter p ~ R~ 3 . On the other hand, if there is a non-trivial matter configuration 
in the bulk so that pd ~ R~ n with n < 3, the cosmological constant, the f(7Z) 
corrections and the effective negative energy density are leading effects. 

The conservation equation (|2.27|) (equivalently (|3.33j )) takes a particularly simple 
form in the low-energy limit in which the term ~ k 2 dominates in the rhs. We obtain 



[p + 3H(p + p)] + e 2 [p d + 3H(p d + p d )] = 0. 



(4.7) 



This the effective four-dimensional description of the energy exchange between the 
brane and the bulk. 



5. Discussion 



The generalisation of the DGP model with a f{TZ) term in the 4D action generates 
novel features in the brane cosmological evolution. An important question concerns 
the presence of regimes of accelerating expansion. The inclusion of higher-derivative 
terms in the gravitational action may lead to such regimes. A typical example is the 



Starobinsky model in four dimensions ||25|| . The evolution equations that we derived 
incorporate the Starobisnky model and its generalizations. The main behaviour 
is more transparent if we neglect the contributions from the matter sectors and 
the spatial curvature. Then, the effective Friedmann equation ( [4.3| ) at low energies 
becomes 



H 2 + 



t2 



kr r 



1 + kr c 



(f-i)H 2 --(nf 



n+HTZf" 



2(l-e 2 )(l + Ar c ) 



(5.1) 
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withft = 6(# + 2# 2 ). 

The branch with €2 = 1 leads to a generalization of the evolution in standard 
f(lZ) gravity, with the higher derivative terms being multiplied by the factor kr c / (1 + 
kr c ). It is easy to check that the evolution equation admits accelerating solutions. 
For / = TZ+TZ 2 / (6E 2 ) we obtain behaviour similar to the Starobinsky model [^3], . 
There are approximate solutions of the form 

aE 2 

H = H l ---{t-t t ) + 0{E 4 ) (5.2) 


TZ = 12H 2 -aE 2 + 0(E 4 ) } (5.3) 

as long as the energy scale E satisfies H 3> E, and a = (1 + kr c )/kr c . For kr c ^> 1, 
we have a = 1 and our setup reproduces exactly the 4D Starobinsky model. In the 
opposite limit kr c < 1, we have a ^> 1 and the regime of accelerated expansion is 
much shorter than in the conventional scenario. 

The branch with €2 = — 1 is the generalization of the self-accelerating branch of 
DGP cosmology || [7], [| 15 1. The acceleration is induced by an effective cosmological 



constant ~ (1 + kr c )/r 2 . If the higher derivative terms are neglected, the constant 
term in the rhs of equation (|5.1|) leads to accelerated expansion with constant H both 
for kr c ^> 1 and kr c <C 1. The second term in the lhs of equation ( |5.1|) generates terms 
proportional to higher powers of H. In general, accelerating solutions are expected 
to exist even when the higher derivative terms are important. In the particular case 
of the Starobinsky model with / = 7Z + 7Z 2 /(6E 2 ), the higher derivative terms give a 
vanishing contribution to the self- accelerating solution. It must be emphasized that, 
within the standard DGP model, the self-accelerating branch has pathologies that 
render its physical signficance questionable flO |. Similar problems may be present in 



the generalized version of the model that we are considering. This issue will be the 
subject of future work. 

Our analysis, which was based on the use of covariant techniques, leads to a 
very intuitive understanding of the effects of matter on the cosmological expansion 
at low energies. The brane matter contributes similarly to matter in conventional 4D 
cosmologies. The effect of bulk matter can be incorporated in a term characterized 
as energy density of the generalized dark radiation. An appropriate pressure can be 
defined for this mirage component, according to equation ( |3.32| ). In terms of these 



quantities the effect of the bulk on the cosmological evolution can be cast in a com- 
pletely 4D form. A conservation equation can be derived (equation Q4.7D ), describing 
the energy exchange between the brane and the bulk, or, in other words, between 
the brane matter and the generalized dark radiation. For the conventional branch 
with 62 = 1, this equation indicates that the total energy of the brane and mirage 
matter sectors is conserved, while the energy density is diluted by the expansion. For 
the self- accelerating branch with e 2 = —1, equation ( |4.7| ) allows the energy of the 
two sectors to grow indefinitely. Another peculiar feature is the appearance of the 
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contribution from the generalized dark radiation with a negative sign in the effective 
Friedmann equation (|4.5|) or (|4.6|) . This feature may lead to a cosmological evolution 



that crosses the phantom divide w = —1 |L5]]. The negative sign is not a consequence 



of the violation of any energy condition in the bulk-brane configuration, as both the 
bulk and brane matter components are assumed to have positive energy density. 

At early times and high energies, the behaviour becomes more complicated be- 
cause of the presence of several unconventional contributions originating in the pecu- 
liarities of RS cosmology, the DGP model and f(lZ) cosmology. Despite the complex- 
ity of the setup, we have provided the equations that describe the general evolution. 
They are the effective Friedmann equation ( |2.25| ) or (LI), and the effective conser- 



vation equation ( p.27| ) or ( |3.33| ). They are applicable to all theories with an induced 
f(7Z) term and an arbitrary number of brane and bulk fields or matter components. 
They provide the basis for the analysis of the features of any specific model. 
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